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We analyse the generation of quantum discord by means of the dynamical Casimir effect in su¬ 
perconducting waveguides modulated by superconducting quantum interferometric devices. We 
show that for realistic experimental parameters, the conditions for the existence of quantum discord 
are less demanding than the previously considered for quantum entanglement or non-classicality. 
These results could facilitate the experimental confirmation of the quantum nature of the dynam¬ 
ical Casimir effect radiation. Moreover, the states with non-zero discord and zero entanglement 
generated by the dynamical Casimir effect are a useful resource for quantum cryptography. 


The Dynamical Casimir Effect (DCE) consists in the 
generation of photons out of the vacuum of a quantum 
field, by means of the modulation of a boundary condi¬ 
tion -e. g. a mirror- at relativistic speeds. The phe¬ 
nomenon was theoretically predicted in 1970 jl]. How¬ 
ever, it was not until 2011 that the first experimental 
observation was reported [2]. The reason for this is that 
accelerating a mirror to relativistic speeds by mechanical 
means is out of experimental reach. The difficulty was 
circumvented by implementing a different kind of mirror. 
In particular, a Superconducting Quantum Interferomet¬ 
ric Device (SQUID) which interrupts a superconducting 
transmission line provides a boundary condition equiv¬ 
alent to a mirror under some experimental conditions. 
Unlike its mechanical counterpart this boundary condi¬ 
tion can be made to move at velocities close to the speed 
of light in the medium, enabling the observation of radi¬ 
ation generated by the DCE. 

A distinctive feature of the DCE radiation is its gen¬ 
uinely quantum nature, due to its origin in the virtual 
particles of the quantum vacuum. In particular, the col¬ 
lected radiation should display entanglement, in stark 
distinction to the one originated from the modulation 
of a thermal source. However, the experimental confir¬ 
mation of the existence of DCE entanglement has proven 
difficult so far. In [3], it is shown that entanglement 
shows up above a critical threshold of thermal noise and 
the conditions for the experimental amplitudes and ve¬ 
locities of the boundary condition necessary to overcome 
that noise are provided. Experimental efforts towards it 
are currently going on in the laboratory. In this work, we 
analyse an alternative strategy to confirm the quantum 
nature of the DCE radiation. We consider the generation 
of a more general form of quantum correlations, such as 
quantum discord. 

Quantum discord HE] has attracted a great deal of 
attention in the last years and represents a paradigm 
shift in the analysis of quantum correlations. In par¬ 
ticular, it has been shown that some quantum states 
that do not possess entanglement still display some gen¬ 
uinely quantum form of correlations -characterised by the 


discord- that are useful for quantum technologies, such 
as remote state preparation d IB], quantum metrology 
0, microwave quantum illumination [TO] and quantum 
cryptography ra¬ 
in this paper, we compute the quantum discord gen¬ 
erated in the experimental setup employed for the obser¬ 
vation of the DCE, namely a superconducting waveguide 
interrupted by a SQUID. Using realistic experimental pa¬ 
rameters, we show that the driving amplitudes and ve¬ 
locities that are required to overcome a certain level of 
thermal noise and generate quantum discord are smaller 
than the ones necessary to generate quantum entangle¬ 
ment and other nonclassicality indicators discussed in [3] . 
In other words, we show that for a given experimental 
value of the driving amplitude and velocity, the critical 
value of temperature above which quantum discord van¬ 
ishes is higher than in the case of quantum entanglement. 
These results could facilitate the experimental confirma¬ 
tion of the quantum nature of DCE radiation. Moreover, 
the kind of states with non-zero discord and zero entan¬ 
glement generated by the DCE in our setup have proven 
to be a useful resource in device-dependent continuous 
variable quantum cryptography [Ill- 

Let us now discuss our model and results in detail. We 
will consider the same experimental setup as in HH- 
The electromagnetic field confined by a superconducting 
waveguide is described by a quantum field associated to 
the flux operator $(x,t)., which obeys the 1+1 D Klein- 
Gordon wave equation, d xx $(x,t) — v~ 2 d t t^{x, t) = 0. 
The field can thus be written in the form 
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where a(co) and b(u>) are the annihilation operators for 
photons with frequency w propagating to the right - 
incoming- and left -outgoing-, respectively. Here we have 
used the notation a(— ui) = a\uj), and k ^ = oj/v is the 
wavenumber, v is the speed of light in the waveguide, and 
Z 0 the characteristic impedance. 
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As shown in mini, for large enough SQUID plasma 
frequency, the SQUID provides the following boundary 
condition to the flux field: 

$(0,i) + L eS (t)d x <f>(x,t)\ x=0 = 0, (2) 

that can be described by an effective length 

L eS {t) = [<S> 0 /2v) 2 /(Ej(t)L 0 ), (3) 

where Lq is the characteristic inductance per unit length 
of the waveguide and Ej(t) = Uj[$ e xt(Q] is the flux- 
dependent effective Joseplison energy. For sinusoidal 
modulation with driving frequency u>d/ 27r and normal¬ 
ized amplitude e, Ej[t ) = Ej[l + e sin u^i], we obtain 
an effective length modulation amplitude (5L e ff = 
where = L e ff(0). If the effective velocity u e ff = 
SL e ffUJd is large enough as compared to v, the emission 
of photons by means of the DCE is sizeable. 

Within this framework, the DCE is analysed us¬ 
ing scattering theory which describes how the time- 
dependent boundary condition mixes the otherwise in¬ 
dependent incoming and outgoing modes [121 ■ I n the 
perturbative regime discussed analytically in minim, 
the resulting output field is correlated at modes with an¬ 
gular frequencies such that w + + = u>d, so 

we can write co± = uid /2 ± Suj, where 6ui is the detuning. 
Introducing the notation a± = a(u>±) and b± = b(ui±), 
the relation between the input and the output fields is 
the following: 

b± = —a± — i —— Uw+oj-at, (4) 

v ' r 

where SL e g^/ui_u) + /v is a small parameter. If we consider 
small detuning, then ~ ~ u>d /2 and 

SL eS ^/uj-UJ + ^ eL eS (0)Lo d _ VeS 

v 2v 2v 

Denoting the small parameter as /, we can write: 

b± = -a± -ifa^. (6) 

Let us consider now the covariance matrix of the sys¬ 
tem V. Using the same convention as in [3] 

Vafi = ^ {RaRfi + RpRa) , 

-which assumes zero displacement- where 

R t = (q-,p-,q+,p+) 


is a vector with the quadratures as elements: 

q± = (b± + tf ± )/y/2 


P± = ~i{b± - b\)/V2. 


Note that the quadratures of the outgoing modes can be 
written in terms of the ingoing modes 

Qo± = (a± +a±)/V2 


and 


Po± = ~i(a± - o±)/V2 
by using Eq. @: 

q± = ~(qo± + Jpot) 

p± = ~(po± + (?) 


We assume that the ingoing modes are in a weakly ther¬ 
mal, quasi-vacuum state characterised by a small fraction 
of thermal photons nV 1 , ?r t i h as is the case for typical GHz 
frequencies and mK temperatures in a superconducting 
scenario. Then the ingoing covariance matrix is 


*>4 


/l + 2n t _ h 0 0 0 

0 1 + 2n t _ h 0 0 

0 0 1 + 2 n‘ + h 0 

\ 0 0 0 l + 2n t + V 


( 8 ) 


Note that since we are considering small detuning 


w+ ~ uj- ~ uJd/2 


and then 


Using Eqs. (J7| and ([8]), we obtain the covariance matrix 
of the outgoing modes 


F = 2 


A C 
C T B 


A = l + 2n t _ h + / 2 (l + 2n£)l, 
B = l + 2n t + h + / 2 (l + 2n t _ h )l, 
C = 2/(1 + n+ + n t })a x . 


(9) 


This is a two-mode squeezed thermal state charac¬ 
terised by the squeezing parameter 2/ and its standard 
form is obtained by just replacing a x by a z in C. The 
main aim of this work is to characterise the quantum cor¬ 
relations of the state described by the covariance matrix 
in Eqs. §• As a measurement of quantum correlations, 
we choose quantum discord which can be exactly com¬ 
puted for two-nrode squeezed thermal states, as recently 
proven by m- In fact, the optimality result of Ref. [13] 
implies that, for these states, the (unrestricted) quantum 
discord coincides with the upper bound of the Gaussian 
discord HEP ] which here takes the form:: 

Cm = h(s/l2)-h(d-)-h(d + )+h( AA_AEIAlh y 

( 10 ) 


and 
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where 


I\ = det A, 1-2 = det B, I 3 = det C , 


and d+, d- are the symplectic eigenvalues, that is 
the eigenvalues of the matrix itt V -where fl = 
(— ia y , 0; 0, —iG y )~ and 

h(x) = {x+] ) ) log 2 (ir + - {x - i) log 2 (z - ^). (11) 

Note that the quantum discord is not symmetric under 
the exchange of modes m , and the expression resulting 
of the exchange of + and — would replace I\ by 12- In 
our case, however, the only effect of the replacement is 
the exchange of n + and n_ , which we are considering to 
be approximately equal in the regime of small detuning. 
Therefore, we can consider just a single expression for 
the quantum discord, Eq. ( |To| . 

Using Eqs. ([9]) and (10), we find that in the perturba¬ 
tive regime: d± = 1/2 — f 2 /2 + n^ 1 and finally 


l'n th 'i 2 

C(V) = max{0, f 2 — —2~}’ (12) 

Therefore, the outgoing state displays nonzero discord as 
long as / > n th /V2- For the sake of comparison with the 
results of [3], we write it explicitly in terms of the driving 
amplitude e, finding that the onset of quantum discord 
occurs at: 


£o 


V2v 

- 1 

Letted 


th 


(13) 


which represents an improvement of \/2 with respect to 
the results reported in [3] for the logarithmic negativity 
mm and other non classicality indicators. 

In Figs. ([!]) and ([i]) we compare the behaviour of VC 
with the logarithmic negativity computed in [3] for realis¬ 
tic experimental parameters. We choose \J~C instead of C 
itself in order to include the same perturbative orders in 
both cases, as explained in [20j. We find that the amount 
of quantum discord is always larger than the entangle¬ 
ment and that the conditions for the existence of nonzero 
correlations are less demanding in the case of discord. In 
particular, for a fixed value of the temperature the dis¬ 
cord is larger than 0 at values of the driving amplitude 
where entanglement is still 0 -Fig. ([!])- or conversely, for 
a fixed value of the driving amplitude the discord is still 
finite at values of the temperature where the entangle¬ 
ment has already vanished -Fig. ©• Therefore, the ex¬ 
perimental conditions required to achieve quantum cor¬ 
relations in the DCE radiation are less demanding than 
those required to achieve entanglement. This could help 
the experimental characterisation of the quantum nature 
of the DCE radiation. Moreover, these results widen the 
applicability of DCE radiation for quantum technologies. 
As a first particular example, the states in Eq. © have 



FIG. 1: Quantum discord \[C -red, dashed- and logarith¬ 
mic negativity -blue, solid- as a function of the normalised 
driving ampitude e. We consider experimental parameters 
v = 1.2 • 10 8 m/s, u>d = 2-7T • 10GHz, L e ff(0) = 0.5mm and 
T = 50mK. Thus the small parameter /< 0.05 is well within 
the perturbative regime, as well as the average numbers of 
thermal photons n th ~ 8 • 10 - ‘ ! . The onset of quantum dis¬ 
cord appears before the entanglement one and the magnitude 
of quantum discord is always larger than entanglement. 



FIG. 2: Quantum discord \[C -red, dashed- and logarithmic 
negativity -blue, solid- as a function of the average number of 
thermal photons n th . We consider experimental parameters 
v = 1.2 • 10 8 m/s, uid = 2-7T • 10GHz, L e ff(0) = 0.5mm and 
t = 0.15. Thus the small parameter / < 0.05 is well within the 
perturbative regime. While quantum entanglement vanishes 
at T ~ 60mK, quantum discord is more robust and survives 
up to T ~ 67mK. 


proven to be the resource in a device-dependent Quan¬ 
tum Key Distribution protocol in the parameter regime 
where they possess zero entanglement and non-zero dis¬ 
cord (see the Supplementary Information of JTj). 

In summary, we have considered the experimental 
scenario of the DCE experiment in a superconducting 
waveguide terminated by a SQUID and analysed the 
quantum correlations of the generated radiation, as char¬ 
acterised by the quantum discord instead of entangle¬ 
ment. We have found that quantum discord is non-zero 
at realistic experimental values of the driving amplitude 
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and temperature where entanglement is zero. Therefore, 
we have extended the experimental parameter regime 
where the quantum nature of the DCE radiation in su¬ 
perconducting scenarios can be assessed. Interestingly, in 
this new parameter regime where discord is different from 
zero while entanglement remains null we have identified a 
first technological application in the realm of continuous 
variable quantum cryptography. A thorough investiga¬ 
tion of the role of DCE quantum discord in quantum 
technologies lies beyond the scope of the current work, 
but promising possibilities might include quantum inter¬ 
ferometric setups (H| or remote state preparation pro¬ 
tocols [7] where quantum discord is a key ingredient. 
This would complement current entanglement-based in¬ 
vestigations on the use of DCE in quantum technologies 
mi- Moreover, our analysis can be extended to differ¬ 
ent quantum platforms. In particular, the question of 
the quantum nature of the phononic DCE radiation in 
Bose-Einstein condensates remains open [22| so a similar 
approach as the one taken here should be of interest in 
that case as well. 
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